Topological Protection of Majorana Qubits 



Meng Cheng,^'^''^ Roman M. Lutchyn/ and S. Das Sarma^ 

^Station Q, Microsoft Research, Santa Barbara, CA 93106-6105 
^Condensed Matter Theory Center, Department of Physics, 
University of Maryland, College Park, Maryland 20748, USA 
* Kavh Institute for Theoretical Physics, University of California, Santa Barbara, CA, 93106 

(Dated: December 19, 2011) 

We study the stability of the topological quantum computation proposals involving Majorana 
fermions against thermal fluctuations. We use a minimal realistic model of a spinless + ipy 
superconductor and consider effect of excited midgap states localized in the vortex core as well as of 
transitions above the bulk superconducting gap on the quasiparticle braiding, interferometry-based 
qubit read-out schemes, and quantum coherence of the topological qubits. We find that thermal 
occupation of the midgap states does not affect adiabatic braiding operations but leads to a reduction 
in the visibility of the interferometry measurements. We also consider quantum decoherence of 
topological qubits at finite temperatures and calculate their decay rate which is associated with the 
change of the fermion parity and, as such, is exponentially suppressed at temperatures well below 
the bulk excitation gap. Our conclusion is that the Majorana-based topological quantum computing 
schemes are indeed protected by the virtue of the quantum non-locality of the stored information 
and the presence of the bulk superconducting gap. 



PACS numbers: 

I. INTRODUCTION 



Topological quantum computation, based on the en- 
coding of quantum information in non-local degrees of 
freedom, provides a promising route to fight quantum 
decoherenc^il^. Due to the presence of environmental in- 
teractions, decoherence is the foremost challenge in any 
conventional quantum computation schemes since quan- 
tum error correction protocols typically have very severe 
constraints on the amount of error that can be corrected 
in a fault-tolerant manner. Topological quantum compu- 
tation (TQC) utilizes topological degeneracy of certain 
low-dimensional systems believed to host non-Abelian 
quasiparticles (anyons): fractional quantum Hall states 
{e.g. at filling factor v = 5/2)^^, certain exotic lattice 
spin system^, and topological superconductor^^^^. The 
latter has received tremendous attention recentljff^HSl 
and a number of possible candidates for topological 
superconductivity has been proposed including stron- 
tium ruthenat el^, t opological insulator-superconductor 
heterostructur^^, semiconductor-superconductor het- 
erostructur^i3t2lla,nd non-centrosymmetric superconduc- 
torP2. In all these systems, non-Abelian Ising anyons are 
realized as Majorana zero-energy modes bound to certain 
to pological d efects and obey non-Abelian braiding statis- 
|-j(,^4l5M23l241_ Given that Majorana zero-energy quasi- 
particles are described by hermitian operators 7 = 7^, 
one can show by constructing a non-local Dirac opera- 
tor c out of two spatially-separated Majorana operators 
7z=i,2 («-e. c = 71 -I- 172) that there are 2"^^ degenerate 
states for a given overall fermion parity with 2n Majo- 
rana zero-energy modes at fixed positions. The infor- 
mation is encoded in the occupied or unoccupied states 
of the non-local Dirac fermion modes. This is a crucial 
concept for the Majorana-based TQC proposals. As long 



as global fermion parity in the system is preserved, one 
can design fault-tolerant quantum computation schemes 
at sufficiently low temperatures. 

In this paper we investigate the effect of finite- 
temperature thermal fluctuations on three key aspects of 
topological quantum computation: quantum coherence 
of the topological qubits, topologically-protected quan- 
tum gates and the read-out of qubits. Since the infor- 
mation is encoded in non-local degrees of freedom of the 
ground state many-body wavefunction, it is important 
to keep the system close to the ground state. How- 
ever, any systems realized in the laboratory are operated 
at a finite temperature T > 0. To prevent uncontrol- 
lable thermal excitations, it is generally accepted that T 
has to be way below the bulk excitation gap. However, 
complications appear when there exist various types of 
single-particle excitations with different magnitudes of 
gaps which can change the occupation of the non-local 
fermionic modes. Note that throughout the paper we 
assume that Majorana fermions are sufficiently far away 
from each other and neglect exponentially small energy 
splitting due to inter-vortex tunneling. The effect of 
these processes on topological quantum computing has 
been discussed elsewher e^ '^^1 Another trivial effect not 
considered in this work is a situation where the fermion 
parity conservation is explicitly broken by the Majorana 
mode being in direct contact with a bath of fermions 
(electrons and holes) where obviously the Majorana will 
decay into the fermion bath, and consequently decohere. 
Such situations arise, for example, in current topological 
insulators where the existence of the bulk carriers (invari- 
ably present due to the unintentional bulk doping) would 
make any surface non-Abelian Majorana mode disappear 
rather rapidly. Another situation that has recently been 
considered in this contexlP3 is the end Majorana mode 
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in a one-dimensional nanowire being in contact with the 
electrons in the non-superconducting part of the semi- 
conductor, leading to a zero-energy Majorana resonance 
rather than a non-Abelian Majorana bound state at zero 
energy. The fact that the direct coupling of Majorana 
modes to an ordinary fcrmionic bath will lead to its de- 
coherence is rather obvious and well-known, and does not 
require a general discussion since such situations must be 
discussed on a case by case basis taking into account the 
details of the experimental systems. In particular, the 
reason the quantum braiding operations in Majorana- 
based systems involves interferometry is to preserve the 
fermion parity conservation. Our theory in the current 
work considers the general question of how thermal fluc- 
tuations at finite temperatures affect the non-Abelian 
and the non-local nature of the Majorana mode. 

We consider a simple model for two-dimensional chi- 
ral Px + ipy superconductor where Majorana zero-energy 
states are hosted by Abrikosov vortices. The quasipar- 
ticle excitations in this system are divided into two cat- 
egories: a) Caroli-de Gennes-Matricon (CdGM) or so- 
called midgap states localized in the vorte x core with en- 
ergies below the bulk superconducting gafpsHH (the gap 
that separates the zero-energy state to the lowest CdGM 
state is called the mini-gap Am)', b) extended states with 
energies above the bulk quasiparticle gap which is de- 
noted by A. The natural question arising in this con- 
text is how these two types of excitations affect topo- 
logical quantum computation using the Majorana zero- 
energy states at finite temperature. This question is 
very relevant in the context of strontium ruthenate as 
well as other weak-coupling BCS superconductors where 
the Fermi energy Ep is much larger than the supercon- 
ducting gap A in which case Am oc A'^ /Ep ^ A. We 
mention in passing here that the semiconductor-based 
Majorana proposal^i^^^^SI ^j^q not have low-lying CdGM 
states because the minigap Am ~ Ap2l due to the small 
Fermi energy in the semiconductor. If the temperature 
is substantially below the minigap, i.e. T <^ Am, obvi- 
ously all excited states can be safely ignored. However, 
such low temperatures with T ^ Am can be hard to 
achieve in the laboratory since for typical superconduc- 
tors A/Ep ^ 10~^ — lO"**. We note that even in the semi- 
conductor sandwich structures the Majorana energetics 
obey the inequality Am < A since in general Ep > A 
even in the semiconductor-based systems in view of the 
fact that typically A ~ 1 K. This makes our considera- 
tion in this paper of relevance also to the semiconductor- 
based topological quantum computing platforms. In this 
paper we investigate the non-trivial intermediate temper- 
ature regime Am ^ 2^ < A. To make this paper more 
pedagogical, we will use a simple physical model that 
captures the relevant physics. We find that the presence 
of the excited midgap states localized in the vortex core 
does not effect braiding operations. However, the midgap 
states do affect the outcome of the interferometry exper- 
iments. 

We also study the quantum dynamical evolution and 



obtain equations of motion for the reduced density ma- 
trix assuming that the finite temperature is set by a 
bosonic bath (e.g. phonons). We find that the qubit 
decay rate A is given by the rate of changing fermion 
parity in the system and is exponentially suppressed (i.e. 
A cx cxp(— A/T)) at low temperatures in a fully-gapped 
Px + "i-Py superconductor. In this context, we make some 
comments about Refs.[?T] claiming to obtain different re- 
sults regarding the effect of thermal fluctuations. 

The paper is organized as follows. In Sec.|lT]we gener- 
alize the notion of non-Abelian braiding to a finite tem- 
perature and show that braiding is not affected by CdGM 
bound states. In Sec.|lTT]we show that the midgap states 
are important for interferometry experiments and gener- 



ally reduce the visibility of the signal. In Sec. IV we study 
the problem of qubit decoherence and effects of thermal 
fluctuations. Finally, we conclude in Sec|V] Some tech- 
nical details are given in the Appendices. 



II. NON-ABELIAN BRAIDING IN THE 
PRESENCE OF MIDGAP STATES 

In this section we address the question of how the 
midgap states affect the non-Abelian statistics at finite 
temperature. The usual formulation of the non-Abelian 
statistics as unitary transformation of the ground states 
does not apply, since at finite temperature the system has 
to be described as a mixed state. We need to generalize 
the notion of the non-Abelian braiding in terms of phys- 
ical observableJ^. This can be done as the following: 
consider a topological qubit made up by four vortices la- 
beled by a = 1, 2,3,4. Each of them carries a Majorana 
zero-energy state, whose corresponding quasiparticle is 
denoted by jao which satisfies 7^q = 1 , 7ao = 7ao ■ There 
are other midgap states in the vortex core which are de- 
noted by dai,i — 1,2, ... ,m. (Actually the number of 
midgap states is huge and the midgap spectrum eventu- 
ally merges with the bulk excitation spectrum. However, 
since we are interested in T <C A, we can choose an en- 
ergy cutoff A such that T ^ A ^ A and only include 
those midgap states that are below A.) It is convenient 
to write dai = 7a,2i-i +*7a,2i, so each vortex core carries 
odd number of Majorana fermions 7ai, i = 0, 1, . . . , 2m. 

Having the notations set up, we now define a gen- 
eralized Majorana operator Fq — i"^ Y[^=o7ai- It is 
straightforward to check that {r^,?;,} = 2Sab- We then 
define the fermion parity shared by a pair of vortices 
'Sab = *ror6- The topological qubit can be uniquely 
specified by a set of measurements of the expectation 
value of the following Pauli matrices & — (ctj., dy^dz)'- 



— S32,i5'j, — I]i3,(T2 — E21- 



(1) 



The non-Abelian braiding can be represented as the 
transformation of (ct). 

Now we list the key assumptions to establish the non- 
Abelian properties of the vortices: 
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1. The fermion parity T,ab is a physical observable that 
can be measured by suitable interferometry exper- 
iments, even at finite temperature. 

2. All the bound states remain localized together with 
the zero-energy state when the vortices are trans- 
ported. Therefore they can be considered as one 
composite system. 

3. The tunneling processes of fcrmions between differ- 
ent vortices and transitions to the gapped contin- 
uum are exponentially suppressed due to the pres- 
ence of the bulk superconducting gap. 

Under these conditions, the only local dynamical pro- 
cesses are the transitions of fermions between the local- 
ized bound states, e.g. scattering by collective excita- 
tions like phonons. However, such processes necessarily 
conserve Fq, therefore also the parities Sad- 

To see this explicitly, the state of the qubit is described 
by the density matrix p(t). Because we are truncating the 
whole Hilbert space to include only those below our cut- 
off A, it is necessary to use the time-dependent instan- 
taneous basi^SH, We can generally consider transitions 
between the various states induced by four-fermion in- 
teractions or coupling to a bosonic bath. To be specific, 
we write down the Hamiltonian of the system: 

H = Ho + ^int. (2) 

Here Hq is the Hamiltonian of the BCS superconductor 
with vortices, whose positions are time-dependent. 
At each moment of time Hq can be diagonalized, yield- 
ing a set of complete eigenbasis which are represented by 
the time-dependent generalization of the aforementioned 
Bogoliubov quasiparticles ■jaoit)i daiit). Hmt describes 
all kinds of perturbations that are allowed under the as- 
sumptions. 

Without going into the details of microscopic calcula- 
tions, we write down the general Lindblad form of the 
master cquatiorP^ governing the time-evolution of the 
density matrix: 

f^^f^- ^[Ht{t),p] + SpS^ - \{S^S, p}. (3) 

The ^ denotes the change of p solely due to the change 

of basis states. Here Ht describes the (effective) uni- 
tary evolution of the density matrix due to transitions 
between different fcrmionic states and the Lindblad su- 
peroperators S corresponds to non-unitary evolution in- 
duced by system-environment coupling. Our assumption 
on the locality of the interactions in the system implies 
that 

[Ht,±ab]^0,[S,±ab]^0- (4) 

The time evolution of the expectation values of cr(i) is 
given by 

d(cr) d _ „ dcr , ^ „ dp 



With it is straightforward to check that 

Ty a[H, p] ^ 0,Ty rT{SpS^ -^{S^S,p}) ^0. (6) 

Therefore we have 

= Tr +Tr .(t) . ^.Tr [.(t)p(t)]. 

(7) 

As we have defined, dt means that all changes come from 
the change in the basis {^ai{t)}- Since after the braiding 
the system returns to its initial configuration, the op- 
erators jia undergo unitary transformations. So if the 
braiding starts ai t — ti and ends at t — tf, we have the 
simple result {a{ti)) — {a{tf)). However, the operators 
r{tf) are different from T{ti). One can easily verify that 
the operators fa satisfy Ivanov's rul^^Ml under braiding 
of vortices a and b: 

fa^ffc,ffc^-f,. (8) 

And the transformation of (ct) is identical to the case 
without any midgap states. In conclusion, in terms of 
physically measurable quantities, the non-Abelian statis- 
tics is well-defined in the presence of excited midgap 
states localized in the vortex core. 

This result is thoroughly non-obvious because it may 
appear on first sight that arbitrary thermal occupancies 
of the mid-gap excited states would completely suppress 
the non-Abelian nature of the system since the Majo- 
rana mode resides entirely at zero energy and not in the 
excited mid-gap states. 

III. INTERFEROMETRY IN THE PRESENCE 
OF MIDGAP STATES 

We now discuss the effect of the midgap states in in- 
terferometry experiments designed for the qubit read- 
ou^ )35tj38| 'pjjgj.g jg number of recent proposals for in- 
terferometry experiments in topological superconduc- 
tord^SMl] Jn this paper, we use an example of the Mach- 
Zehnder interferometer proposed by Grosfeld and Ster 
based on Aharonov-Casher (AC) effect. In this proposal, 
a Josephson vortex (fluxon) is driven by supercurrent 
to circumvent a superconducting island. The fluxon ap- 
pearing at the interface between two topological -\- ipy 
superconductors carries a zero-energy Majorana modes, 
and behaves as a non-Abelian anyon. Therefore, the vor- 
tex current around the central superconductor is sensitive 
to the topological content of the enclosed superfluid. (We 
refer the reader to Ref. \42^ for more details.) Indeed, 
the vortex current is proportional to the total tunneling 
amplitude: 

Ja « \{tLUL+tLUR)\'fo)\^ 

= KlP + liflP + 2Re{tltR{'^o\UL'UR\'^o)} (9) 
= KlI' + l^fll' + 2Re{tltfle''^*^(*o|M|*o)}. 



4 




FIG. 1: Mach-Zehnder interferometer proposed in Ref. [35] 
for topological qubit detection. Due to the Aharonov-Casher 
effect, the vortex current is sensitive to the charge enclosed. 
Long Josephson junction between two topological supercon- 
ductors carries allows for Josephson vortices (fluxons) that 
carry Majorana zero-energy modes. 



Here |^o) is the initial state of the system and Ul and 
Ur are the unitary evolution operators for the fluxon 
taking the two respective paths, ^pac is the Aharonov- 
Casher phase accumulated by the fluxon: i^ac — T^Q/e. 
Here Q is the total charge enclosed by the trajectory of 
the fluxon, including the offset charge Qoxt set by ex- 
ternal gate and the fermion parity rip of the low-energy 
fermionic states: 



Q = Qcxt + eup 



(10) 



M encodes the transformation solely due to the braid- 
ing statistics of the non-Abelian fluxon around n non- 
Abelian vortices. 

If the superconducting island contains no vortices, then 
M = 1 and the interference term is solely determined 
by the AC phase. The magnitude of the vortex current 
shows an oscillation: 



Jv — JvQ 



1 + C COS 



ttQ 

e 



(11) 



Here ( is the visibility of the interference. 

When n is odd, there is no interference because 
Mj^'o) and |4'o) have different fermion parity, implying 
(^'o|M|^o) = 0. The vortex current becomes indepen- 
dent of the charge encircled. Therefore, the disappear- 
ance of the interference can be used as a signature of the 
non-Abelian statistics of the vortices. 

We now consider a situation where the non-Abelian 
fluxon has midgap states other than the Majorana bound 
state. The internal state of the fluxon then also depends 
on the occupation of these midgap states. As we have 
argued in the previous section, as far as braiding is con- 
cerned the non-Abelian character is not affected at all by 
the presence of midgap states. So the interference still 
vanishes when there are odd numbers of non-Abelian vor- 
tices in the island. On the other hand, when there are 
no vortices in the island, transitions to the midgap states 
can significantly reduce the visibility of the interference 
term C. 



To understand quantitatively how the visibility of the 
interference pattern is affected by the midgap state, let us 
consider the following model of the fluxon. Since we are 
interested in the effect of midgap states, we assume there 
is only one midgap state and model the probe vortex by 
a two-level system, or spin 1/2, with the Hilbert space 
{|0),|1)}. Here |1) denotes the state with the midgap 
state occupied. We also assume that the charge enclosed 
by the interference trajectory Q = so we can neglect 
the AC phase. The Hamiltonian is then given by 



H = \L){L\(g)HL + \R){R\(E)He 



(12) 



where Hl,r is given by 



k k 

(13) 

Here rj = L,R. Notice that we assume the bath couples 
to the fluxon locally so we introduce two independent 
baths for L and R paths. The unitary evolution at time 
t is then factorizable: 



Uit) = \L){L\ (g, UL{t) + \R){R\ ® Unit). 



(14) 



Given initial state /5(0) = Ppnth ^ Ps ^ Pbath, we can 
find the off-diagonal component of the final state p{t) = 
If (t) p{0)U^ (t) , corresponding to the interference, as 



Tr 



Pbath 



= Tr PsTri[pbath,L?^L(i)]Trij[pbath,flC^i?(^)] 



(15) 



Now we evaluate Wrj{t) = Tr^[pbath,jjC^i)(i)] (notice W'^ 
is still an operator in the spin Hilbert space). We drop 
the 77 index in this calculation. First we switch to in- 
teraction picture and the evolution operator U{t) can be 
represented formally as U{t) — Texp{— i J^^ dt' Hi{t')} 
where 

H,(t) = ^5fc(a+e*^*/2+a-e-*^*/2)(4e^"''*-fafce-'"^*). 

k 

Following the derivation of the master equation for the 
density matrix, we can derive a "master equation" for 
W{t) under the Born-Markovian approximation: 



dW 
"dT 



= -7(n+ 1/2 + 0- J2)l^, 



(17) 



where7 = irj^k 9kH^k- ^),n = 7 ^ X]fc 5fc'^fe'5('^fc-^)■ 
Therefore, the visibility of the interference, propor- 
tional to the trace of W, is given by 



C oc Tr[W{t)ps] cx e 



-■ynt —--ynL/v 



(18) 



We notice that the model we have used is of course a 
simplification of the real fluxon. We only focus on the 
decoherence due to the midgap states and assume that 
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only one such state is present. In reality, there could be 
many midgap states which lead to a stronger suppression 
of visibility. 

The above interferometor is able to detect the existence 
of non-Abelian vortices which requires that the Joseph- 
son vortex (i.e. fluxon) also has Majorana midgap states. 
To fully read out a topological qubit, one needs to mea- 
sure the fermion parity of the qubit. This can also be 
done using interferometry experiments with flux qubits, 
essentially making use of the AC effect of Josephson vor- 

ticcs43 44. 

Another relevant question is whether the thermal ex- 
citations of the (non-Majorana) midgap states localized 
in the vortex core have any effects on the interferome- 
try. Since the interferometry is based on AC effect where 
vortex acquires a geometric phase after circling around 
some charges, one might naively expect that the inter- 
ferometric current might depend on the occupation of 
the midgap states due to the charge associated with the 
midgap states (i.e. Qn — e J dr(|u„p — |z)„p)). The 
situation is more subtle, however, once one takes into 
account the screening effect due to the superfluid con- 
densate. The kinetics of the screening process is beyond 
the scope of this paper. However, assuming equilibrium 
situation more careful inspection discussed in the Ap- 
pendix [B] shows the geometric phases acquired by the 
Josephson vortices only depends on the total fermion par- 
ity in the low-energy midgap states (even if they are not 
Majorana zero-energy modes) and the offset charge set 
by the external gate voltage. 



then reads 



IV. DEPOLARIZATION OF QUBITS AT FINITE 
TEMPERATURE 



We now study the coherence of the topological qubit 
itself. From our discussion on the effect of bound states 
in the vortex core, it is clear that decoherence only oc- 
curs when the qubit is interacting with a macroscopically 
large number of fermionic degrees of freedom, a fermionic 
bath. An example of such a bath is provided by the con- 
tinuum of the gapped quasiparticles in the topological 
superconductor. Once the Majorana fermion is coupled 
to the bath via a tunneling Hamiltonian, the fermion 
occupation in the qubit can leak into the environment, 
resulting in the depolarization of the qubit. It is then 
crucial to have a fully gapped quasiparticle spectrum to 
ensure that such decoherence is exponentially small, as 
will be shown below. 

To study the decay of a Majorana zero mode, we con- 
sider two such modes, 71 and 72, forming an ordinary 
fermion c = 71 + ij2 ■ The gapped fermions are coupled 
locally to 7i, without any loss of generality. The cou- 
pling is mediated by a bosonic bath. The Hamiltonian 



k I 

-"i^gkiliidi + dk){d\ + di). 

kl 



(19) 



Here is the annihilation operator of the gapped 
fermions with quantum number k and energy Sk- di is 
the annihilation operators of the bosonic bath. 

The density matrix of the system evolves according 
to the equation of motion p = —i[H,p\. Since we are 
interested in the qubit only, we will derive the master 
equation for the reduced density matrix pr, tracing out 
the bosonic bath and the gapped fermions. 



dt 



= -A[p,-7i(-l)"/),(-l)"7i] 



(20) 



where 



'M + niinl + ime.-LoO. (21) 



kl 



Here n{ = l/(e^'=/^ + l),n'l = l/(e"'/'^ - 1) are the 
Fermi and Bose distribution functions. The derivation of 
( 20 1 is presented in the Appendix |Aj Notice that at low 
temperatures T <C A, due to energy conservation, both 



and 



n{ are suppressed by the Gibbs factor e 



-A/T 



Therefore, the rate A ^ e^^^"^ . 

Then the polarization of the qubit (cr^) — Tr[azP,:] sat- 
isfies dt{cz) = —'^^{o'z)- Therefore the lifetime of the 
topological qubit is given by Ti ^ A^^. Physically, this is 
reasonable since we introduce tunneling term between the 
Majorana fermion and the gapped fermionic environment 
so the fermion parity of the qubit is no longer conserved. 
It is expected that A is determined by the exponential 
factor e^^l^ when T <^ A. Therefore, this provides a 
quantitative calibration of the protection of the topolog- 
ical qubit at finite temperature. In the high-temperature 
limit T ^ A, the distribution function scales linearly 
with T so the decay rate is proportional to T. This is 
quite expected since T ^ A, the gap does not play a 
role. We note that a recent work by Goldstein and Cha- 
morPil studying the decay rate of Majorana zero modes 
coupled to classical noise essentially corresponds to the 
high-temperature limit of our calculation T ^ A and, as 
such, does not apply to any realistic system where the 
temperature is assumed to be low, i.e. T <ti A. In fact, 
in the trivial limit of A <C T, the Majorana decoherence 
is large and weakly temperature dependent because the 
fermion parity is no longer preserved and the fermions 
can simply leak into the fermionic batlP^. By defini- 
tion, this classical limit of T ^ A is of no interest for 
the topological quantum computation schemes since the 
topological superconductivity itself (or for that matter, 
any kind of superconductivity) will be completely absent 
in this regime. Our result makes sense from the qual- 
itative considerations: quantum information is encoded 
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in non-local fermionic modes and changing fermion par- 
ity requires having large thermal fluctuations or external 
noise sources with finite spectral weight at frequencies 
a; ~ A. Furthermore, it is important to notice that 
such relaxation can only occur when the qubit is cou- 
pled to a continuum of fermionic states which renders 
the fermion parity of the qubit undefined. Intuitively, 
the fermion staying in the qubit can tunnel to the con- 
tinuum irreversibly, which is accounted for by the proce- 
dure of "tracing out the bath" in our derivation of the 
master equation. It is instructive to compare this result 
to a different scenario, where the zero-energy fermionic 
state is coupled to a fermionic state (or a finite number 
of them) instead of a continuum. In that case, due to 
hybridization between the states the fermion number os- 
cillates between the two levels with a period (recurrence 
time) determined by the energy difference IS.E between 
them. The expectation value of the fermion number (or 
spectral weight) in the zero-energy state is depleted and 
oscillatory in time, but will not decay to zero. 

The above derivation can be straightforwardly gener- 
alized to > 2 Majorana fermions, each coupled locally 
to gapped fermions and bosonic bath. 

^ = - ^ A, [p, - ^,(-l)"-p,(-l)".^,] , (22) 
i=i 

The depolarization of the qubit can be calculated in the 
same fashion. 



V. CONCLUSION AND DISCUSSION 

We study quantum coherence of the Majorana-based 
topological qubits. We analyze the non-Abelian braid- 
ing in the presence of midgap states, and demonstrate 
that when formulating in terms of the physical observ- 
able (fermion parity of the qubit), the braiding statistics 
is insensitive to the thermal occupation of the midgap 
states. We also clarify here the conditions for such topo- 
logical protection to hold. Our conclusion applies to the 
case of localized midgap states in the vortex core which 
are transported along with the Majorana zero states dur- 
ing the braiding operations. If there are spurious (e.g. 
impurity-inducecP^Hll) midgap bound states spatially lo- 
cated near the Majorana zero-energy states but are not 
transported together with them, they could strongly af- 
fect braiding operations. For example, during braiding 
the fermion in the qubit has some probability (roughly 
determined by the non-adiabaticity of the braiding op- 
eration) to hybridize with the other bound states near 
its path leading to an error. If the disorder is weak and 
short-ranged, such low-energy states are unlikely to oc- 
cur unless the bulk superconducting gap is significantly 
suppressed at some spatial points (e.g. vortices) as it is 
well-known that for a single short-range impurity the en- 
erg y of s uch a bound state is close to the bulk excitation 
gap^SM' Thus, well-separated impurity-induced bound 



states are typically close to the gap edge and would not 
affect braiding operations. If the concentration of impuri- 
ties is increased, then it is meaningful to discuss the prob- 
ability distribution of the lowest excited bound state in 
the system'^^. The distribution of the first excited states 
determing the minigaps depends on many microscopic 
details (e.g. system size, concentration of the disorder). 
Since the magnitude of the minigaps is system-specific, 
one should evaluate the minigap for a given sample. As 
a general guiding principle, it is important to reduce the 
effect of the disorder which limits the speed of braiding 
operations. However, we note here that physically mov- 
ing anyons for braiding operations might not be necessary 
and there are alternative measurement-only approaches 
to topological quantum computatior!^ where the issue of 
the low- lying localized bound states is not relevant. 

In this paper we also consider the read-out of topolog- 
ical qubits via interferometry experiments. We study the 
Mach-Zehndcr interferometer based on Aharonov-Casher 
effect and show that the main effect of midgap states 
in the Josephson vortices is the reduction of the visibil- 
ity of the read-out signal. We also consider the effect of 
thermal excitations involving midgap states of Abrikosov 
vortices localized in the bulk on the interferometry and 
find that such processes do not effect the signal provided 
the system reaches equilibrium fast enough compared to 
the tunneling time of the Josephson vortices. 

Finally, we address the issue of the quantum coher- 
ence of the topological qubit itself coupled to a gapped 
fermionic bath via quantum fiuctuations. We derive the 
master equation governing the time evolution of the re- 
duced density matrix of the topological qubit using a sim- 
ple physical model Hamiltonian. The decoherence rate 
of the qubit is exponentially suppressed at low tempera- 
tures T <C A. Since topological protection assumes that 
fermion parity in the superconductor is preserved, our 
result is very intuitive. 

We conclude that the Majorana-based qubits are in- 
deed topologically well-protected at low temperatures as 
long as the experimental temperature regime is well be- 
low the superconducting gap energy. 
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Appendix A: Derivation of the Master Equation 

In this appendix we derive the master equation for the 
reduced density matrix. The density matrix of the whole 
system evolves according to the equation of motion: 



dp 
dt 



= t[Hi,p]. 



(Al) 



Notice that we will be working in interaction picture in 
the following. Here the coupling Hamiltonian 



where 



Hi{t) ^ i^gknfikit)<j)i{t). 

kl 



(A2) 



(A3) 



Assume the coupling between the qubit and the bath 
is weak, we integrate the equation of motion for a time 
interval At: 



At 



1 

At 



dti / dt2TrB[^/(ti),[i?/(t2),/5(t2)]]. 

(A4) 

The first-order term vanishes due to the fact that {(t>{t)) = 
{fl{t)) — 0. Now we make the Born approximation for the 
bath: assume that the bath is so large that it relaxes very 
quickly to thermal equilibrium. The density matrix of 
the whole system can be factorized as p{t) — Pr{t) ® Pb- 
Here the bath includes with the gapped fermionic bath 
and the bosonic bath. 



The commutator on the right-hand side of ( A4 1 can be 
evaluated: 



TrB[i?/(ti)Ji//(t2),/5(t2)]]« [/5rW-7(-l)"A-W(-l)"7]{(%(ti)%(i2))(0^ 

(A5) 

I 



The factor (—1)" appears because of the anti- 
commutation relation between fermionic operators. The 
correlators of the bath are easily calculated: 

{m{h)f,k{t2)) = 4e--(*i-*^) + (1 - n{)e-^^^(*i-*=) 

{Mti)Ut2)) = nl'e*"'^*^-*^' + {n\ + l)e--H*i-*^) 

(A6) 

Performing the integral over ti and t2, we finally arrive 

at 

^ = -A[p, - 7i(-l)"A-W(-l)"7i] . (A7) 
Here A is given by 
\ = 2Y,gl, [(1 - + ni{n\ + - ui). (A8) 

kl 



Appendix B: Geometric Phases Generated by 
Midgap Fermions 

In this Appendix we derive the geometric phase gener- 
ated by the midgap fermions, relevant to the interferom- 
etry experiments involving Josephson vortices (fluxons). 
We follow here the formalism developed in the context of 
AC effect for flux qubits— ■. 

We assume that a superconducting island with several 
midgap fermionic states, labeled by is coupled to a 
flux qubit. In the low-energy regime well below the bulk 
superconducting gap and the plasma frequency, the only 



degrees of freedom of this system are the superconduct- 
ing phase (f) and the midgap fermions. We also assume 
that the phase varies slowly so the fermionic part of the 
system follows the BCS mean-fleld Hamiltonian with su- 
perconducting phase 4>. 

We want to know the geometric phase associated with 
vortex tunneling in the presence of midgap fermions. It 
can be derived by calculating the transition amplitude 
Afi associated with a time-depedent phase (p = 0(t): 

Af,^{cl,f\QfU{tf,U)Q\\cj,,), (Bl) 

where \(f)) denotes the BCS ground state with supercon- 
ducting phase (j) and 4>f—(j)i — 2wtt. = nm(^m)"" 
note the occupation of the midgap fermionic states with 
rim = 0,1. 

The midgap fermionic operators dj„ are explicity ex- 
pressed in terms of Bogoliubov wavefunctions Um and 

dlit) = e— dr [7.,„(r)7/;t(r)e^^/2+«„(r)^(r)e-^^/2]. 

(B2) 

Therefore, 

U{tf,U)dl{U)uHtf,U) = dl(t/)e— (B3) 
So the transition amplitude is evaluated as 

= {Pf\QfU{tf,U)Qlu\tf,U)U{tf,UM,) 

(B4) 
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We conclude that the geometric phase is precisely irwri = 
^{4>f — (f>i) Physically this reflects the fact that one 
fermion is "half" of a Cooper pair. The vortex tunneling 
causes the phase of the Cooper pair condensate changes 
by 27r and correspondingly the fermionic states obtain tt 



phases. Notice that the phase '^„i£m{tj — t-i) is simply 
the overall dynamical phase of the whole system due to 
its finite energy and does not contribute to the interfer- 
ence at all. 



^ A. Y. Kitacv, Ann. Phys.(N.Y.) 303, 2 (2003). 

^ C. Nayak, S. H. Simon, A. Stern, M. Freedman, and 

S. Das Sarma, Rev. Mod. Pliys. 80, 1083 (2008). 
^ P. Bonderson, S. Das Sarma, M. Preedman, and C. Nayak, 

arXiv: 1003.2826 (2010). 
" G. Moore and N. Read, Nucl. Pliys. B 360, 362 (1991). 
5 C. Nayak and F. Wilczek, Nucl. Phys. B 479, 529 (1996). 
® S. Das Sarma, M. Freedman, and C. Nayak, Phys. Rev. 

Lett. 94, 166802 (2005). 
^ N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 
* A. Kitaev, Physics- Uspekhi 44, 131 (2001). 
^ D. A. Ivanov, Phys. Rev. Lett. 86, 268 (2001). 
1° F. Wilczek, Nat. Phys. 5, 614 (2009). 
" A. Stern, Nature 464, 187 (2010). 
^2 M. Franz, Physics 3, 24 (2010). 
13 B. G. Levi, Physics Today 64, 20 (2011). 
1* S. Das Sarma, C. Nayak, and S. Tewari, Phys. Rev. B 73, 
220502 (2006). 

1^ L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 (2008). 
1® A. Cook and M. Franz, Phys. Rev. B 84, 201105 (2011). 
^"^ J. D. Sau, R. M. Lutchyn, S. Tewari, and S. Das Sarma, 

Phys. Rev. Lett. 104, 040502 (2010). 
1® J. Alicea, Phys. Rev. B 81, 125318 (2010). 
1^ R. M. Lutchyn, J. D. Sau, and S. Das Sarma, Phys. Rev. 

Lett. 105, 077001 (2010). 
^° Y. Oreg, G. Refael, and F. von Oppen, Phys. Rev. Lett. 

105, 177002 (2010). 

21 X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Phys. Rev. B 82, 
184516 (2010). 

22 M. Sato and S. Fujimoto, Phys. Rev. B 79, 094504 (2009). 

23 P. Bonderson, V. Gurarie, and C. Nayak, Phys. Rev. B 83, 
075303 (2011). 

2* J. Alicea, Y. Oreg, G. Refael, F. von Oppen, and M. P. A. 

Fisher, Nat. Phys. 7, 412 (2011). 
2^ M. Cheng, R. M. Lutchyn, V. Galitski, and S. Das Sarma, 

Phys. Rev. Lett. 103, 107001 (2009). 
2^ M. Cheng, R. M. Lutchyn, V. GaHtski, and S. Das Sarma, 

Phys. Rev. B 82, 094504 (2010). 

J. D. Sau, C. H. Lin, H.-Y. Hui, and S. Das Sarma, 
arXiv:1103.2770 (2011). 
2® C. Caroli, P. de Gennes, and J. Matricon, Phys. Lett. 9, 
307 (1964). 

2^5 N. B. Kopnin and M. M. Salomaa, Phys. Rev. B 44, 9667 
(1991). 



3° J. D. Sau, R. M. Lutchyn, S. Tewari, and S. Das Sarma, 
Phys. Rev. B 82, 094522 (2010). 

31 G. Goldstein and C. Chamon, arXiv: 1107.0288 (2011). 

32 A. R. Akhmerov, Phys. Rev. B 82, 020509 (2010). 

33 M. Cheng, V. Galitski, and S. Das Sarma, Phys. Rev. B 
84, 104529 (2011) 



34 
35 



G. Lindblad, Commun. Math. Phys. 48, 119 (1976). 
S. Das Sarma, M. Freedman, and C. Nayak, Phys. Rev. 
Lett. 94, 166802 (2005). 
3® A. Stern and B. I. Halperin, Phys. Rev. Lett. 96, 016802 
(2006). 

3'' P. Bonderson, A. Kitaev, and K. Shtengel, Phys. Rev. Lett. 

96, 016803 (2006). 
3® P. Bonderson, K. Shtengel, and J. Slingorland, Ann. Phys. 

(N. Y.) 323, 2709 (2008), ISSN 0003-4916. 
3^ A. R. Akhmerov, J. Nilsson, and C. W. J. Beenakker, Phys. 

Rev. Lett. 102, 216404 (2009). 
1° J. D. Sau, S. Tewari, and S. Das Sarma, Phys. Rev. B 84, 

085109 (2011). 

E. Grosfeld, B. Seradjch, and S. Vishveshwara, Phys. Rev. 
B 83, 104513 (2011). 
■12 E. Grosfeld and A. Stern, Proc. Natl. Acad. Sci. USA 108, 
11810 (2011). 

"3 J. D. Sau, S. Tewari, and S. Das Sarma, Phys. Rev. A 82, 

052322 (2010). 

•1* F. Hassler, A. R. Akhmerov, C.-Y. Hon, and C. W. J. 
Beenakker, New J. Phys. 12, 125002 (2010). 
D. J. Clarke and K. Shtengel, New J. Phys. 13, 055005 
(2011). 

R. M. Lutchyn, T. D. Stanescu, and S. Das Sarma, 
arXiv:1011.5643 (2011). 

R. M. Lutchyn, T. D. Stanescu, and S. Das Sarma, Phys. 
Rev. Lett. 106, 127001 (2011). 
*® T. D. Stanescu, R. M. Lutchyn, and S. Das Sarma, Phys. 
Rev. B 84, 144522 (2011). 

A. V. Balatsky L Vekhter, and J.-X. Zhu, Rev. Mod. Phys. 

78, 373 (2006). 

^" Q.-H. Wang and Z. D. Wang, Phys. Rev. B 69, 092502 
(2004). 

^1 P. W. Brouwer, M. Duckheim, A. Romito, and F. von Op- 
pen, Phys. Rev. Lett. 107, 196804 (2011). 

*2 p. Bonderson, M. Freedman, and C. Nayak, Phys. Rev. 
Lett. 101, 010501 (2008). 



